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Coherent feedback control of quantum systems has demonstrable advantages over measurement- 
based control, but so far there has been little work done on eoherent estimators and more specifically 
eoherent observers. Coherent observers are input the coherent output of a specified quantum plant, 
and are designed such that some subset of the observer and plant’s expectation values converge in 
the asymptotic limit. We previously developed a class of mean tracking (MT) observers for open 
harmonic oscillators that only converged in mean position and momentum; Here we develop a class 
of covariance matrix tracking (CMT) coherent observers that track both the mean and covariance 
matrix of a quantum plant. We derive necessary and sufficient conditions for the existence of a 
CMT observer, and find there are more restrictions on a CMT observer than there are on a MT 
observer. We give examples where we demonstrate how to design a CMT observer and show it 
can be used to track properties like the entanglement of a plant. As the CMT observer provides 
more quantum information than a MT observer, we expect it will have greater application in future 
coherent feedback schemes mediated by coherent observers. Investigation of coherent quantum 
estimators and observers is important in the ongoing discussion of quantum measurement; As they 
provide estimation of a system’s quantum state without explicit use of the measurement postulate 
in their derivation. 


I. INTRODUCTION 

Quantum engineering has seen rapid growth in the 
last two decades. Physicists, mathematicians and en¬ 
gineers have been working in unison to control a num¬ 
ber of diverse systems in the quantum regime [ll-[8|. 
Quantum control involving feedback has become partic¬ 
ularly topical [ikl, as using information gained from 
a system can lead to more stable operation of a control 
protocol Quantum feedback can be split into 

two paradigms: measurement-based and coherent feed¬ 
back. Measurement-based feedback involves some mea¬ 
surement step in the feedback loop ini, unfortunately 
measurements of quantum systems are typically slow and 
noisy as they involve coupling small quantum systems to 
macroscopic read out devices. Coherent feedback on the 
other hand is feedback where the controller and system 
are coupled directly without a measurement step [IT 
[ 3 . The advantage is that the time scales of the con¬ 
troller and system can be made very similar as they are 
on the same scale. But beyond this practical advantage, 
there is increasing evidence that retaining the coherence 
of the feedback signal provides an intrinsic advantage 
over measurement-based feedback 111 Elis 113 . 

Coherent feedback is still a relatively new paradigm, 
and as such it lacks many of the tools commonly used 
in classical, or for that matter, other quantum feedback 
schemes. In particular, there are still only a limited num¬ 
ber of options for coherently estimating a state within a 
feedback loop. It is well established classically that es¬ 
timation using Kalman filters can provide improved per- 
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formance over direct feedback schemes [HI, and similar 
demonstrations have been performed for measurement- 
based quantum feedback [H|- Unfortunately, traditional 
techniques do not appear to be applicable to coherent 
feedback due to difficulties with quantum conditioning 
onto non-commutative subspaces of signals. Instead, the 
closest option is so-called coherent quantum observers 
[23|. We previously developed a class of coherent quan¬ 
tum observers (as shown in Fig. [T]), which can estimate 
the observables of linear and bilinear quantum plants de¬ 
scribed by quantum stochastic differential equations (QS- 
DEs) in the sense of mean values, independent of any 
additional quantum noise in the observer [nil. We 
have proved that MT coherent observers can always be 
found, consistent with the laws of quantum mechanics, if 
the plant is detectable. In some cases the estimation of 
mean-values is sufficient and feedback can be improved 
with a coherent quantum observer. However, in many 
cases, the energy, correlations and indeed entanglement 
of the observed system may be the target of control or 
needed for feedback, for which the coherent quantum ob¬ 
server would not provide a reliable estimate. To remedy 
this issue, we propose to develop a modified coherent ob¬ 
server to track mean values, variances and correlations, 
namely, the CMT coherent observer. 

In general, a CMT coherent observer outperforms a 
MT coherent observer in several respects. For instance, 
a CMT coherent observer allows us to achieve the most 
similar quantum state to that of the plant. Furthermore, 
it is well established that for a two-mode linear Gaussian 
system, the quantum correlations is completely charac¬ 
terized by the first and second moments [HI, and thus 
entanglement can be mimicked by the utility of a CMT 
coherent quantum observer in this situation. Therefore, 
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one can conclude a CMT observer can provide a better 
estimate in most cases. Nonetheless, we find that the er¬ 
ror convergence rate of a CMT coherent observer can not 
be made arbitrarily high, plus we cannot guarantee that 
CMT coherent observers exist for systems where mean 
value coherent observers exist. 
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FIG. 1. A quantum plant and the corresponding coherent 
quantum observer in a cascade arrangement 

The paper is organized as follows. We begin in Section 
[Uby presenting the linear quantum state space model for 
open harmonic oscillators in the Heisenberg picture. In 
Section nni we briefly discuss quantum plants and (MT) 
coherent quantum observers. In Section lEl we analyze 
the existence of CMT coherent observers, and show the¬ 
orems which tell us how to construct CMT observers to 
be consistent with the laws of quantum mechanics. This 
is followed by numerical simulations in Section O which 
illustrate the design and performance of CMT observers. 
Section IVTI provides some concluding remarks and future 
research directions. The mathematical notation we use 
is defined in Appendix [Al 


II. OPEN HARMONIC OSCILLATORS AND 
LINEAR QSDES 

The dynamics of an open quantum system are u niquely 
determined by the parametrization (S^L^H) [26l-[^. 
The self-adjoint operator H is the Hamiltonian describ¬ 
ing the self-energy of the system. The unitary matrix S is 
a scattering matrix, and the column vector L with opera¬ 
tor entries is a coupling vector. S and L together specify 
the interface between the system and the fields. In the 
physics literature, it is common practice to describe open 
quantum systems using a master equation for a density 
operator p, and it can easily be obtained from the triple 
indeed, we have 

dp = {i [p, H] + C* {p)) dt (1) 

where C* (p) = L'^pL^ — ^L^L'^p — ^pL^L'^ (notation 
defined in Appendix EJ and we assume natural units are 
being used. Given an operator X defined on the initial 
Hilbert space H, its Heisenberg evolution is defined by 

dX = {C {X) - i [X, H]) dt + dW^S^ [X, L] 

+ [L^,X]SdW + tr[{S^XS-X)dA^]. (2) 

With 

C{X) = ^L^X,L] + ^[L^,X]L, (3) 


which is called the Lindblad superoperator (Note >C* (•) 
is the adjoint superoperator of >C(-)). The operators W 
are defined on a particular Hilbert space called a Fock 
space F. When the fields (the number of fields is n^) 
are in the vacuum states, these are the quantum Wiener 
processes which satisfy the Ito rule 

dWdW^ = In^dt. 

Input field quadratures W + W'^ and —i (W — W**) are 
each equivalent to classical Wiener processes, but do not 
commute. A field quadrature can be measured using ho¬ 
modyne detection [^, . The gauge processes A^^, are 

input signals to the system as well. 

We assume there is no interaction between different 
fields, and thus hereafter we assume S to be the identity 
matrix without loss of generality [T5|. This assumption 
eliminates the first time on the right hand side of Eq. ©■ 
To be specific, 

dX ={C{X) -i[X,H])dt 

+ ^{[X,L]-[X,L^])dWi 

-^-{[X,L] + [X,L^)dW2, (4) 


with 



-i{W-W^) 


The quadrature form of the output fields is given by 


' dYi ' 


_dY2 _ 



-i{L-L^) 


dt T 


dWi 

dW2 


( 5 ) 


In this work we focus on open harmonic oscilla¬ 
tors. The dynamics of each oscillator is described by 
two Hermitian operators position qi and momentum 
Pi, which satisfy the canonical commutation relations 
[Qi^Pj] = 2i6ij where Sij is the Kronecker delta. For 
our purposes, it is convenient to collect the position 
and momentum operators of the oscillators into an n- 
dimensional column vector x(t), defined by x {t) = 

(Qi (t) ,Pi (t), ?2 (t) ,p2it),...,qn (t) ,Pn ■ In this 
case the commutation relations can be re-written as: 


X (t) X (t)"^ — (^X (t) X (t)^^ = 2iGr, 


(6) 


where 0^ = J with J = 

Harmonic oscillators, in particular, are defined by hav¬ 
ing a quadratic Hamiltonian of the form H = ^x^Rx 
with R being a symmetric matrix, and a coupling 

operator of the form L = Ax with A being a 
matrix (here n, and Uy are positive even numbers). 
A special property of open harmonic oscillators is that 
the differential equations governing x{t) are linear. If we 
use an n^-dimensional column vector y (t) to incorporate 
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all the quadratures of the output fields then, based on 
Eqs. dH) and ([5]), the dynamics of a set of open harmonic 
oscillators can be described by the following linear QS- 

DEs [l3: 


We assume the quantum plant is some system of open 
harmonic oscillators with a set of Ap^Bp^Cp and Dp ma¬ 
trices which are known but we are unable to change. The 
linear QSDEs (see Section El) for the plant is then: 


dx (t) = Ax (t) dt + Bdw (t), (7a) 

dy (t) = Cx (t) dt + Ddw (t) (7b) 


where ^4, 5, C, D are ^ ]^riyXn 

matrices respectively defined in terms of H and L as 

follows: 


A = 2Gn (i? + S(AtA)), 
B = 2i0„[-At A^]r„, 


c = 

D = 


[ Uy 


T^y 

2 

1 - 

el 

X 

O 

A + A# 

0 ^ riw 

2^2 

Tny 

2 

—iA + iA'^ 


X Tly^ j 5 


with 



Em = Pml^ ^ 



1 i 

1 —i 


(8a) 

(8b) 

(8c) 

(8d) 


and the symbol {tu is a positive even number) de¬ 
notes an m X m permutation matrix defined so that if we 
consider a column vector a = [ ai a 2 • • • Um] , then 

Pm^ — [ ^3 ' ' ' 1 ^2 ^4 ' ' ' ] 

In this work we are primarily interested in engineering 
the A^B^C and D matrices rather than deriving them 
from H and L. When engineering, instead of using Eqs. 
dH]) we instead typically use the so-called physical realiz¬ 
ability conditions: 

AQji ~ 0, (9a) 

BD^ = enC^eny, (9b) 

These are algebraic constraints, independent of H and L, 
which the coefficient matrices A^B^C and D must obey 
for them to correspond to a physically realizable quan¬ 
tum system. They were originally derived by requiring 
the canonical commutation relations of x{t) {y{t)) must 
hold for all times, a property enjoyed by open physical 
systems undergoing an overall unitary evolution l29| . 
But since it has been proven: given a set of A, B^C and 
D matrices that satisfy Eqs. m a corresponding H and 
L can always be found that satisfy Eqs. (jHl) (e.g. see 
Theorem 3.4 in 0). 


III. QUANTUM PLANTS AND COHERENT 
QUANTUM OBSERVERS 

The primary goal of this work is to create a eoherent 
quantum observer which asymptotically tracks the ob¬ 
servables of some arbitrary quantum plant [Ullllll. 


dxp (t) = ApXp (t) dt + Bpdwp (t), (10a) 

dyp {t) = CpXp (t) dt + Dpdwp (t) (10b) 


where Ap, Bp^ Cp are and ma¬ 
trices respectively (here and Uy^ are positive even 


numbers), and Dp = 


yp '^yp^{'^wp '^yp^ 


Eurther- 


more, A^, Bp^ Cp and Dp satisfy the following physical 
realizability conditions 


Ap 0 n^ + P Pp^riwpPp — O 5 (11a) 

B,D^ = en^c^eny^. (11b) 

As shown in Eig. [U we take the quantum output signal 
of the plant and directly fed it into the coherent quantum 
observer [^ . [ 3 ^ . 

An (MT) eoherent quantum observer is another system 
of quantum harmonic oscillators which we engineer such 
that the system variables track those of the quantum 
plant asymptotically in the sense of mean values. The 
coherent quantum observer is driven by the output of the 
quantum plant directly; No measurement is involved. A 
coherent quantum observer has equations of the form 


dxo (t) = {Ap - KCp) Xo (t) dt Kdyp (t) + Bodwo (t), 

(12a) 

dyo (t) = CoXo {t) dt + Do[ dyp (t)’^ dwo {t)^ ] ^ (12b) 


where: the na,-dimensional column vector Xq (t) de¬ 
notes the “estimate” of Xp{t); A, Bo are MA^^^yp ^ 
]^na;xn«;^ matrices respectively; and Do is given by Do = 


^ny„x{nyj,+n„^-ny„) ■ Note that the System de- 

scribed by Eqs. m must also satisfy the following phys¬ 
ical realizability conditions 


{Ap - KCp) + 0„, (Ap - KCp)^ 

+ KQnyK^ + BSny,^ = 0, (13a) 

[K Bo]Dl = Q^^ClQny^ (13b) 

which put restrictions on K and B^ [2^. In the case of 
7 ^ 0, the algebraic constraints Eqs. m indicate that 
an additional quantum noise signal vuQ{t) is needed. 

We use Pp (t) and po {t) fo denote the first moments of 
the plant and the observer respectively, i.e.. 


Pp (t) = {Xp {t)) , 
Po {t) = {Xo {t)) ■ 


The equations of motion for the first moments of the 
plant and the observer are: 


Ap (^) — ^ph'p (^) 5 (14a) 

Po {t) = {Ap - KCp) Po (t) + KCpPp {t ). (14b) 
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Now we define (t) = (t) — fio {t) as the error which 

gives the difference between the first moments of the 
plant and the corresponding observer. And according 
to Eqs. (na) , it evolves as 

{t) = (Ap - KCp) Bp, (t). (15) 

converges to zero asymptotically if and only if Ap — 
KCp is Hurwitz [ 2 ^. Hurwitz here means that all the 
eigenvalues of Ap —KCp have strictly negative real parts, 
and hence lim eCt) = 0 . 

t^oo 

Thus given a quantum plant described by Eqs. m, 
the coefficient matrices of a MT coherent quantum ob¬ 
server described by Eqs. (nni) are designed such that 

1 . {Ap — KCp) is Hurwitz; 

2 . The system described by Eqs. m corresponds to 
an open quantum harmonic oscillator. 

Eurt her more, a MT coherent quantum observer can 
always be found with arbitrary rates of error convergence 
(proportional to the real parts of eigenvalues of Ap — 
KCp) for a detectable plan t The term detectable 

comes from classical control [^, , and it means that all 

modes of the plant are either observable or stable. Where 
observability means only given the outputs the state of 
a mode can be determined in finite time. Whether a 
plant is detectable or not can be .j udg ed entirely from 
the Ap and Cp matrices (e.g., see [^, for details). 
Detectability then forms a sufficient (but not necessary) 
condition for a MT coherent observer to exist. 

A MT coherent observer is limited in that it only tracks 
the mean values of the plant. The covariance matrix is 
not guaranteed to match between the plant and the MT 
coherent observer. This means that important proper¬ 
ties, e.g. the entanglement or energy of the plant, may 
not be correctly estimated. We aim to remove these limi¬ 
tations and create a CMT coherent observer whose quan¬ 
tum state matches the plant completely in the asymptotic 
limit. 


IV. CMT COHERENT OBSERVERS FOR OPEN 
HARMONIC OSCILLATORS 

As an extension of a MT observer, here we create a 
CMT observer which also tracks the covariance matrix 
of a quantum plant. We require that a CMT observer 
also tracks the mean values, thus every CMT observer is 
also a MT observer (but not vice versa). 

Let us use Ep {t) and (t) to denote the covariance 
matrix of the plant and the observer respectively, and 


Epo denotes the cross variance: 

(t) =\ {^P (t) (t) + {xp (t) xl (tC) 

- {Xp (t)){xl (t)), 

(t) {xo (t) xl it) + {Xo (t) x'^ 

-{xAmxUt)), 

Spo (t) = {Xp (t) xl {t)) - {Xp {t)){x'l {t)). 

The evolutions for the correlation matrices of the plant 


and observer are given by: 

(^) —^p'^p (t) + '^p (t) + ^p^p ^ (16a) 

^po {t) =^p^po {t) + {t) {^p ~ KCp) 

+ Ep (t) {KCp)^ + BpK^, (16b) 

^po it) =Ko it) + {Ap - KCp) {t) 

+ KCp^p it) + KBl, (16c) 

So {t) = {Ap - KCp) So {t) + So {t) {Ap - KCp)^ 

+ KCp^po {t) + Sj, {t) {KCp)^ 

+ KK^ + BoBl, (16d) 


where Tip{t), Eo (t) and T^po{t) are real matrices with 
Ep (t) and Eo (t) nonnegative. The difference between 
Ep (t) and E^ (t) is {t) = Ep (t) — E^ (t). Then a CMT 
coherent quantum observer is defined as 

Definition 1 Given a system described by Eqs. m , « 
system described by Eqs. / f71|) is a CMT coherent quan¬ 
tum observer for the system described by Eqs. m tf 

1. The system described by Eqs. m is a MT coher¬ 
ent quantum observer for the system described by 
Eqs. (03); 

2. The covariance matrix of the observer described 
by Eqs. / f7H) tracks that of the plant described by 
Eqs. nU} asymptotically, i.e., 

lim Ep (t) — Eo (t) = lim {t) = 0 . 

t^oo t^oo 

Our main theorem which concerns the existence of a 
CMT coherent quantum observer of the form (HU) is pre¬ 
sented below. 

Theorem 1 There exists a CMT coherent quantum ob¬ 
server described by Eqs. m for a quantum plant de¬ 
scribed by Eqs. m if and only if 

1. Ap — KCp is Hurwitz; 

2. The following identity 

lim {Eo 0 Eo - Ep ^ Ep) 

X (s/ 4^2 — In2x ^ A — A ^n2x ) 

X vec {BB^) = 0 


(17) 
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holds. Here 


Ep — [ Ifix ] 5 

Eo = [ On^ Irix ] 5 


whieh gives Eq. m 

Finally, Eqs. m must eorrespond to an open har- 
monie oseillator, whieh requires that the physieal realiz¬ 
ability condition given by Eqs. /fl^) should hold fldi . [H/. 


and the coefficient matrices of a joint plant- 
observer system are given by 


A = 
B = 


Ap 0 

KCp Ap - KCp \ ’ 

Bp 0 

KDp Bo \ ’ 


3. The system described by Eqs. m is physically re¬ 
alizable. 


Proof 1 First of all, in order to ensure the convergence 
of Cp {t), Ap — KCp must be Hurwitz. 

The covariance matrix for the joint plant-observer sys- 


^po (t) 
So (t) 

lowing Lyapunov differential equation 


tern denoted E (t) = 




if) 

po (t) 


satisfies the fol- 


t (t) = AT, (t) + T (t) + BB'^. (18) 

Note that 

Ep (t) = EpY] (t) Ep , 

So {t)=EoT{t)Ej 

and thus 

vec (Up (t)) = Ep0 Epvec (H (t)), (19) 

vec {Eo (t)) = Eo^ EoYec {E (t)). (20) 


We have found necessary and sufficient conditions for 
the existence of a CMT observer. But it is still a challeng¬ 
ing task to construct a CMT observer by solving Eqs. (HZD 
and m- Thus we consider a special case where it is eas¬ 
ier to construct a CMT observer. Specifically, we assume 
Ap is Hurwitz. Any plant with a unique steady state has 
an Ap matrix which is Hurwitz. 

The primary advantage of Ap being Hurwitz, is that 
we can guarantee the existence of steady states values for 
all the covariance matrices (i.e., lim Ep{t) = 0). Solving 

t^oo 

Eqs. (p!6|) in steady state gives: 

{Ap - KCp) es + es {Ap - KCp)^ 

+ KCp {Tp - Tpo) + {Tp - Ejj {KCp)'^ 

+ BpB^ - KK^ - BoB^ = 0 (22) 

where the steady state = lim {t). 

t^oo 

Eurthermore when Ap — KCp is Hurwitz (as required 
for a CMT observer), then = 0. Substituting = 0 
to Eq. ([22]) gives 

KCp {Tp - Tpo) + {Tp - Sjj {KCp)^ 

+ BpB^ - KK^ - BoBl = 0 (23) 

in steady state. 


By using the Laplace transform £ (•) to Eq. [TM) . we can 
obtain 


C (vec {T {t)))= {shnl - A - A® ) 


-1 


vec 


{BBn 


-h vec(i; (0)) 


Theorem 2 Assume the quantum plant described by 
Eqs. nUi) is detectable with Ap Hurwitz. The system de¬ 
scribed by Eqs. m is a CMT coherent quantum observer 
for the plant described by Eqs. m if and only if 

1. Ap — KCp is Hurwitz; 


then 


2. The following matrix inequality 


E (vec (So {i) - Tp {t))) = 

{Eq ® Eq Ep ® Ep) ^ 5 / 41^2 ^n 2 x ^ ^ ^0 ^n 2 x ) 


vec 


{BB'^ 


+ vec(S(0)) 


KCp {Tp - Tpq) + {Tp - Epo)^ {KCp)^ 

+ BpBl -KK^ -iKQr^^^K^ 

-i{A- KCp) - iQnx {A - KCp)'^ h 0 (24) 


Since a CMT observer has the property that lim (t) = 

t^oo 

0, we have to require that all the poles of C {vec {ej: (t))) 
are located on the left side of the s-plane. Or equivalently, 

lim5£ (vec (e^ (t))) =hm {Eo 0 Eo — Ep 0 Ep) 

s^O s^O 

X (s/4„2 - ®A- A®In 2 x)~^ 

X vec {BB^) 

( 21 ) 


holds, where {Ep — Epo) is the unique solution to 
the following Sylvester equation 

Ap {Ep — Epo) + {Ep — Epo) {Ap — KCp) 

+ HpHj - BpK^ = 0. (25) 

Assuming the two conditions above hold, the coupling op¬ 
erator characterizing the interaction between the observer 
and additional boson fields is then given by Lo = 


=0 
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where is any x Ux eomplex matrix sueh that 
AlAo = - - KCp) -^iA- KCpf On. 

~ ~ ^po) ©n^; 

- J©n. (Sp-Epof {KCpfen. 

+ (26) 

Proof 2 Smce plant deserihed by Eqs. m is de- 
teetable, one ean always find K to make Ap — KCp Hur- 
witz. With the assumption of Ap being Hurwitz, and ae- 
eording to Eq. /f^) . Bq must satisfy 


based on Eq. 

Therefore, Aq is any x Ux eomplex matrix sueh 
that 

AlA, = - ^0„^ (A - KCp) -fiA- KCp)^ Qn. 

~ ~ ^po) ©n^ 

-^©n. (Sp-Spo)^(KCp)^0„, 

- \Qn.BpBlen. + \en.KK^en. 

^ 0 

and viee versa. Eq. [2^ ean then be derived using the 
identity -©n^©n^ = /n^. 


B,Bl =KCp (Sp - Epo) + (Sp - Sjo) (KCp)^ 

+ BpBl - KK^, (27) 

and the eorresponding physieally realizability eondition is 

=-{Ap- KCp) 0„, - 0„, {Ap - KCp)^ 
-KOn^^K^ (28) 

Therefore Bq ean be determined based on Eqs. and 

In aeeordanee with the physieal form of an open har- 
monie oseillator deseribed by Eqs. m with Lo = AqXo, 
Bo is given by lEMM]) 

Bo = 2i©n^ [ — AJ AJ ] (29) 

Here is defined in SeetionUA 

By using the form of Bo given in Eq. we ean 

obtain that 

BoBl = -40„, SR (aJA o) 0„,, (30) 

then 

SR (AtA,) = - \en.KCp {Ep - Epo) 0„, 

{Ep-Ep,)^{KCp)^en. 

-\e„^BpB^en. + \en^KK^en^, ( 31 ) 

due to Eq. 

Similarly, we have 

~ (^ 0 ^ 0 ) ©na;5 (32) 

then 

9 (At Ao) = - ^0„^ (A - KCp) -fiA- KCp)^ 0„, 

+ ^0„,770„,^77^0„, (33) 


As studied in a MT coherent quantum observer 
can always be found if the plant described by Eqs. m is 
detectable. However, as we intend to track the covariance 
matrix of a linear quantum plant using coherent observers 
at the same time, not all values of K that make Ap — KCp 
Hurwitz are applicable to the design of a CMT coherent 
observer. Indeed, there are systems where mean value 
coherent observers exist but CMT observers can not be 
constructed. It is worth mentioning that Bo can be 0 
if no additional noise is needed to ensure the physical 
realizability of an observer described by Eqs. m ■ 


V. APPLICATIONS AND EXAMPLES 

In this section, we present some numerical examples to 
illustrate the design and performance of CMT coherent 
quantum observers. We also compare the behavior of an 
MT vs. CMT observer. 


A. CMT observers vs. MT observers for a 
single-mode quantum harmonic oscillator 


In this example we consider tracking a single-mode 
Gaussian system. Consider an optical parametric oscil¬ 
lator as the linear quantum plant given by 

_ r -0.4 0 

axp- ^ 0 -0.6 

dyp = Xpdt -\- dwp (34b) 


Xpdt — dwp (34a) 


where Ap = 


/2. 


-0.4 0 

0 - 0.6 


, Bp = —I 2 and Cp = Dp = 


If we choose K to be 3 / 2 , then using Eq. (EHl) one 
1 0 


can choose Bo = 
quantum observer 


0 -2 


to construct a MT coherent 
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However, in this case, according to Eq. (EH) we have 


BoB 


T 

o 


-1.6842 0 

0 -2.2857 


(35) 


which is negative, and therefore a CMT coherent observer 
cannot be designed with K = 3 / 2 . 

Alternatively, one can set K = I 2 . First, we can cal¬ 


culate the steady state 


1.1111 0 
0 0.9091 


using Eq. (ESI. Then by substituting K and E^ — E^o 
to Eq. dll]), we find the Eq. holds. Applying the 
Cholesky decomposition, one can determine 


Ao = 


0.6742 0.7416i 
0 0.0745 


It is thus that 


Bo =2i02 [ -aJ ] r4 


-1.4832 0 0 0.1491 

0 -1.3484 0 0 


Also, we choose the initial covariance matrix for the joint 
plant-observer system as 


S(0) = 


1.1/2 0 
0 2/2 


which corresponds to a Gaussian separable joint state 
pBj . The initial amplitudes are (0) = [ 1 1 ] and 

/io(0)= [0 0]^. 

We can calculate E^ (t) and Eo (t) explicitly by using 
the Laplace transform, and 


es (i) = Sp (t) - So (i) 


2 7 14 f 

-h " 




We investigate the convergence of the covariance ma¬ 
trices between the plants and the coherent observers by 
plotting the Frobenius matrix norm of the covariance er¬ 
ror matrix \\ej:{t)\\F = ^Tr[e|] against time in Fig. O 
When ||eE(t)||F = 0 we can be certain eE(t) = 0, and 
hence the covariance matrices of the plant and observer 
are identical. We can see the CMT observer is perform¬ 
ing as expected. The matrix E^ {t) is tracking E^ (t) 
asymptotically as time goes to infinity. 

We also investigate the quantum correlations between 
the plant and the CMT observer. For Gaussian two¬ 
mode systems, entanglement is completely quantified by 
the smallest symplectic eigenvalue (t) of the partially 
transposed state, and the joint state is entangled if and 
only if Z 2 _ (t) < 1 [ 2 ^, [25|. In Fig. [2] we plot the smaller 
symplectic eigenvalue as a function of time. We find that 
the plant and the CMT observer eventually become en¬ 
tangled as depicted by the dash-dot line in Fig. [2l 



FIG. 2. Plot of the behavior of a CMT observer with the 
solid line and the dash-dot line corresponding to ||eE(t)||F 
and 17- (t) (see Eq. (13) in [ 2 ^) respectively. The joint system 
is initialized in a Gaussian separable state. 


As the CMT observer tracks both the first and second 
moments of the plant, and the quantum state is Gaus¬ 
sian, we expect the quantum state of the CMT observer 
to be identical to that of the plant. This is not guaran¬ 
teed to be the case for the MT observer that only tracks 
the means. We compare the performance of the CMT 
and MT observer in this regard by plotting the quantum 
fidelity between the observer and the plant as a function 
of time in Fig. [3l Quantum fidelity is widely used to 
quantify how close two mixed states are m m. For 
Gaussian states, the fidelity between two states can be 
calculated analytically (see Eq. (7) in [33|). In this pa¬ 
per, we use F (t) to denote the fidelity, and the closer 
F (t) is to 1 the more similar the two sates are to each 
other. In Fig. [S] the state of a CMT observer (with 


Bo 


-1.4832 0 0 0.1491 

0 -1.3484 0 0 


) gets closer to the 


plant state compared to a MT coherent observer (with 


Bo 


1 0 
0 2 


), as anticipated. 


B. Entanglement tracking of a two-mode quantum 
harmonic oscillator using a CMT observer 


In this example we consider a linear quantum plant 
which consists of two oscillators that are initially separa¬ 
ble but eventually become entangled. The initial covari¬ 
ance matrix for the plant is E^ (0) = ^ 2 / 


and 
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FIG. 3. Plot of the fidelity F (t) (see Eq. (7) in Q) as a func¬ 
tion of time with the solid line and the dash-dot line corre¬ 
sponding to a CMT observer and a MT observer respectively. 


its evolution is governed by the linear QSDEs 


dxp = 


d'Up — 


-0.4 

0 


0 

0 

-0.1 

6 

0 

1 

0 


Ip 

0 

1 


0 

- -1 

0 

0 

0" 

0 

-1 

0 

0 

1 

0 

1 

0 

0 

1 

0 

2_ 

1 0 

-1 

0 - 


0 1 

0 - 

-1 

X 

0 0 

-2 

0 

0 0 

0 - 

-1 _ 



0 

0 

0 

- 1.6 


dw> 


pi 


Xpdt 


Xpdt + dw. 


ip. 


(36a) 


(36b) 


vZ. {t) converges to (t) asymptotically, as expected. 
This confirms that even quantum correlations inside the 
two-mode Gaussian plant can be tracked by the CMT 
observer. This allows for control of the plant based on 
quantum characteristics that were unavailable with a MT 
observer. 



EIG. 4. Plot of an entanglement measure (smallest symplectic 
eigenvalue) of the observer (t) (the solid line) and the plant 
(t) (the dash-dot line) as a function of time. 


C. Failed tracking of the covariance matrix of a 
singe-mode quantum harmonic oscillator 


Consider a plant with the following linear QSDEs 


We now design a CMT coherent observer for the plant 
described by Eqs. p6|). One can choose the observer gain 
K to be 


K = 


"0.2 

0 

-0.1 

0 ■ 


0 

0.05 

0 

-0.1 

(37) 

0.6 

0 

-0.1 

0 

0 

0.4 

0 

-0.1 _ 



Then we find that Eq. ([24]) holds, and thus a CMT ob¬ 
server can be constructed according to Theorem [2] with 
(Ao is not unique) 

" 0.5167 0.5952i -0.2914 -0.1887r 
0 0.0571 -0.0167i 0.1343 

0 0 0.9316 0.4887i ' 

0 0 0 0.027 

We initialize the observer to Eo (0) = 2/4, which is dif¬ 
ferent to the plant initial condition, but still separable. 

We now confirm that the entanglement between the 
oscillators of the plant is correctly tracked by the CMT 
observer. In Fig. |Hwe plot the smallest symplectic eigen¬ 
value of the partially transposed state of both the ob¬ 
server (t) and the plant (t) as a function of time. 


dXp = 


^Vp — 


-1 1 
1 -1 

V2 0 
0 V2 


Xpdt + 



Xpdt + dwp 


0 

-^/2 


dwp, 


(38a) 

(38b) 


For the quantum plant given by Eqs. (|38|) (Note that Ap 
is not Hurwitz), no matter what values we choose K and 
Bo to be, Eq. (HB cannot be satisfied. It is thus that a 
CMT coherent quantum observer can never be designed 
for the plant given by Eqs. (|38|) . 

There do exist plants which cannot be tracked by a 
MT observer. For instance, certain undetectable plants 
cannot be tracked even in the mean values sense. How¬ 
ever, in this example, a MT observer can be constructed 
even though a CMT observer cannot be. Specifically, we 
can choose K = I 2 to make Ap — KCp Hurwitz, and 
then Bo is determined as shown in [^. Therefore, in 
this case, we are only able to approach the mean val¬ 
ues of the quantum plant without tracking its covariance 
matrix. This demonstrates that there are additional con¬ 
straints when constructing a CMT observer compared to 
a MT observer. 
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VI. DISCUSSION AND OUTLOOK 

We have created a CMT observer that tracks both the 
mean values and covariances of a system of linear quan¬ 
tum oscillators. We first emphasize that the CMT ob¬ 
server (and the previously developed MT observer) does 
not just have the same steady state as the system it is 
tracking, but it is also correlated to the plant’s state. 
This is shown in Fig. 2 for the CMT observer, where we 
see it is entangled with the plant. 

In future work it will be investigated whether an ob¬ 
server can be used in a feedback loop to better control 
the behavior of a plant than would otherwise be possible 
with direct feedback. It has already been shown that co¬ 
herent feedback has advantages over measurement-based 
feedback. We now expect that coherent observer medi¬ 
ated feedback will have an advantage over direct feedback 
in the same sense a classical observer mediated feedback 
has an advantage over direct feedback. 

Furthermore, we expect that the CMT observer will be 
more useful than a MT observer as it tracks the covari¬ 
ance matrix of the plant as well. Important properties of 
a quantum system, such as its energy, entanglement and 
other quantum correlations are a function of the system’s 
covariance matrix rather than merely its mean values. 

Our work on a CMT coherent observer also raises some 
interesting fundamental questions with regard to engi¬ 
neering quantum systems in comparison to classical sys¬ 
tem. Classifying what plants can or cannot be tracked 
with a MT observer appears to be identical to classical 
observer theory. Namely a K must be found such that 
Ap — KCp is Hurwitz. There is well established clas¬ 
sical theory which then relates this requirement to no¬ 
tions such as observability and detectability [3l|, [32|, . 

A CMT observer, on the other hand, has additional 
requirements which are fundamentally quantum in ori¬ 
gin. Namely, Eq. dSll) must be satisfied in addition to 
Ap — KCp being Hurwitz. It remains an open question 
on how to interpret this additional requirement and if the 
classical notions of observability and detectability can be 
appropriately extended when discussing the tracking of a 
quantum plant’s covariance matrix. As we are attempt¬ 
ing to copy the entire quantum state of the plant with 
a CMT observer (unlike a MT observer), there may be 
some connection between these additional requirements 
and the no-cloning theorem (sH, • 

Outside of quantum engineering, the design and im¬ 
plementation of a CMT observer also looks to provide 
some insight into quantum measurement. When the out¬ 
put of the plant is measured, an optimal estimate of 
the quantum state of the plant can be calculated using 
the Belavkin-Kalman filter (also referred to as stochas¬ 
tic trajectories) [9[. However, research suggests the sit¬ 


uation becomes much more complicated when there is 
no measurement step. It has been proven the Belavkin- 
Kalman filter fails in the presence of a fully quantum 
non-commutative output signal and furthermore 

measurement-based Kalman filters are challenging to be 
realized efficiently with quantum hardware [^. The 
CMT observer is the first coherent method of providing 
an estimate of the full quantum state of a plant. Note we 
never invoked the measurement postulate when deriving 
the CMT observer. It is entirely derived in the framework 
of open systems. Creating and better understanding es¬ 
timators for quantum systems which do not explicitly 
require the measurement postulate is an important part 
of further refining our understanding of quantum mea¬ 
surement. 
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Appendix A: Notation 

In this paper the asterisk is used to indicate the Hilbert 
space adjoint of an operator A, as well as the complex 
conjugate = x — iy oi d, complex number z = x iy 
(here, i = and x^y are real). Real and imaginary 

parts are denoted by 5R (z) = and ^ (z) = 

respectively. The conjugate transpose of a matrix 
A = {aij} is defined by A'^ = {ci'ji}^ The conjugate 
= {a^j} and the transpose A^ = {aji} of a ma¬ 
trix is defined so that A^ = (A^)*^ = (A^)^. det (A) 
denotes the determinant of a matrix A, and tr (A) rep¬ 
resents the trace of A. vec(A) denotes the vectoriza- 
tion of a matrix A. ||A||^ denotes the Erobenius norm, 
i.e., IIA11^ = y^tr^AbA). The mean value (quantum 
expectation) of an operator X in the state p is de¬ 
noted by (A) = Ep [A] = tr(pA). The commutator 
of two operators A, Y is defined by [A, Y] = XY — YX. 
The anticommutator of two operators A, Y is defined by 
{A, K} = XY -j- YX. The tensor product of operators 
A, Y defined on Hilbert spaces H, G is denoted X 0 Y, 
and is defined on the tensor product Hilbert space HgG. 
In (n ^ N) denotes the n dimensional identity matrix. On 
(n G N) denotes the n dimensional zero matrix. 
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